NAVAL  POSTGRADUATE  SCHOOL 
MONTEREY,  CALIFORNIA 


THESIS 


Thesis 
Y8226 


EFFECTS  OF  SHALLOW  WATER  BOTTOM  INTERFACE 

ROUGHNESS  AND  VOLUME  FLUCTUATIONS  ON 

BROADBAND  PULSE  RESOLUTION 


by 


Mei-Chun  Yuan 


June,  1996 


Principal  Advisor: 
Second  Reader: 


Kevin  B.  Smith 
Anthony  A.  Atchley 


Approved  for  public  release;  distribution  is  unlimited. 


DUr  LIBRARY 

NAV  RADUATE  SCHOOL 

MONTEREY  GA  93943-5101 


REPORT  DOCUMENTATION  PAGE 


Form  Approved  OMB  No  0704-0188 


Public  reporting  burden  lor  this  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instruction,  searching  existing  data  sources, 
gathering  and  maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this 
collection  ol  information,  including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports.  1215  Jefferson 
Davis  Highway,  Suite  1204,  Arlington.  VA  22202-1302.  and  to  the  Office  of  Management  and  Budget,  Paperwork  Reduction  Project  (0704-0188)  Washington  DC  20503 


AGENCY  USE  ONLY  (Leave  blank) 


REPORT  DATE 
June  1996 


REPORT  TYPE  AND  DATES  COVERED 
Master's  Thesis 


TITLE  AND  SUBTITLE  EFFECTS  OF  SHALLOW  WATER  BOTTOM  INTERFACE 
ROUGHNESS  AND  VOLUME  FLUCTUATIONS  ON  BROADBAND  PULSE 
RESOLUTION 


AUTHOR(S)  Mei-Chun  Yuan 


FUNDING  NUMBERS 


PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESSi  I  S  , 
Naval  Postgraduate  School 
Montcrcx.  C'A  93943-5000 


PERFORMING 
ORGANIZATION 

Kl ■!'<  )l<  I  NUMBER 


SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESSi  ES> 


SP(  )NS(  )RJN(  i/M(  )NIT(  )RIN(  i 
AGENCY  REPORT  NUMBER 


SUPPLEMENTARY  NOTES  The  views  expressed  in  this  thesis  are  those  of  the  author  and  do  not  reflect  the  official 
policy  or  position  of  the  Department  of  Defense  or  the  U.S.  Government. 


12a     DISTRIBUTION/A VATLABLLrTY  STATEMENT 

Approved  for  public  release;  distribution  is  unlimited. 


12b. 


DISTRIBUTION  CODE 


1 3       ABSTRACT  (maximum  200  words) 

Typical  acoustic  propagation  in  shallow  water  environments  is  dominated  by  bottom-interacting 
paths  The  effects  of  rough  bottom  interfaces  and  sediment  volume  fluctuations  are  investigated  using 
model  simulations  A  numerical  study  of  low-frequency  (-200  Hz)  broadband  pulse  propagation  is 
presented  and  several  characterizations  of  bottom  factors  are  examined  In  particular,  the  variations  of 
the  interface  rms  roughness  and  a  volume  fluctuation  strength  constant  on  the  time  resolution  of  the 
broadband  pulse  are  investigated.  It  is  shown  that  interface  roughness  is  the  dominant  factor  in  the 
degradation  of  pulse  resolution. 


14.     SUBJECT  TERMS  Underwater  Acoustic  Propagation,  PE  Model,  Beamforrrnng,  Pulse  Resolution. 
Scattering. 


15.     NUMBER 
PAGES     67 


OF 


16.     PRICE  CODE 


17.      SECURITY  CLASSIFICA- 
TION OF  REPORT 
Unclassified 


SECURITY  CLASSIFI- 

CATION OF  THIS  PAGE 
Unclassified 


19.     SECURITY  CLASSIFICA- 
TION OF  ABSTRACT 
Unclassified 


20.     LIMITATION 
ABSTRACT 
UL 


OF 


NSN  7540-01. 280-5500 


Standard  Form  298  (Rev.  2-89) 

Prescribed  by  ANSI  Std.  239-18  298-102 


11 


Approved  for  public  release;  distribution  is  unlimited. 

EFFECTS  OF  SHALLOW  WATER  BOTTOM  INTERFACE  ROUGHNESS  AND 
VOLUME  FLUCTUATIONS  ON  BROADBAND  PULSE  RESOLUTION 


Mei-Chun  Yuan 

Lieutenant  Commander,  Republic  of  China  Navy 

B.S.,  Chung  Cheng  Institute  of  Technology  -  1987 

Submitted  in  partial  fulfillment 
of  the  requirements  for  the  degree  of 

MASTER  OF  SCIENCE 

IN 

ENGINEERING  ACOUSTICS 

from  the 

NAVAL  POSTGRADUATE  SCHOOL 

June  1996 


IDLEY  KNOX  LIBRARY 

TGRADUATE  SCHOOL 
,    EREY  CA  93943-5101 


ABSTRACT 
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I.  INTRODUCTION 

Over  the  past  few  years,  the  interest  in  littoral,  shallow  water  environments  has 
increased  significantly  The  success  of  many  naval  systems  in  these  environments  depends 
on  our  ability  to  predict  and  interpret  shallow  water  acoustic  propagation  The  sound 
speed  profile  in  shallow  water  typically  consists  of  an  upper  mixed  layer  below  which  the 
sound  speed  decreases  to  the  bottom  Typical  shallow-water  environments  are  found  on 
the  continental  shelf  with  approximate  depths  of  200  m  The  wave  propagation  is 
downward  refracting,  that  is,  bottom-interacting  paths  dominate  in  long-range  propagation 
as  depicted  in  Figures  1  and  Figure  2  Etter  (1991)  summarized  the  effect  of  the  sea  tloor 
as 

•  forward  scattering  and  reflection  loss  (but  complicated  by  refraction  in  the 
bottom), 

•  interference  and  frequency  effects; 

•  attenuation  by  sediments; 

•  noise  generation  at  lower  frequencies  due  to  seismic  activity, 

•  backscattering  and  bottom  reverberation 

The   purpose   of  this   thesis   is   to   quantify   the   influence   of  bottom   variability   on   the 
resolution  of  pulse  transmissions  in  simple,  shallow  water  environments 

Rough  bottom  interfaces  can  be  modeled  by  assuming  some  realistic  form  of  the 
spectral  shape.  The  forward  scatter  of  the  acoustic  field  can  then  be  treated 
deterministically  Fluctuations  of  the  sediment  properties,  such  as  porosity,  density, 
permeability,  and  shear  strength,  determine  the  sound  scattering  from  the  bottom  volume. 
Yamamoto  (1995)  summarized  the  relationship  between  spectral  strength  and  sediment 
type  We  used  this  information  to  generate  random  two-dimensional  wave-number 
spectra  of  the  velocity  variability  in  order  to  simulate  a  realistic  ocean  sediment  volume. 
These  techniques  are  described  in  more  detail  in  Chapter  II. 
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Typical  sound  speed  profile  in  shallow  water 
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Figure  1.  Typical  sound  speed  profile  in  shallow  water. 


Ray  trace  in  shallow  water 
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Figure  2.  Ray  trace  in  shallow  water. 


Wave  propagation  in  an  inhomogeneous  medium  is  defined  by  the  wave  equation. 
The  solution  can  be  obtained  by  several  approximate  methods,  for  example,  Wentzel- 
Kramers-Brillouin  (WKB)  approximation,  ray  approximation  or  parabolic  equation 
approximation  The  last  method  was  used  throughout  this  thesis  study.  Over  the  past  few 
decades,  ocean  acoustic  models  have  been  developed  and  improved  to  analyze  measured 
data  and  simulate  sound  propagation  in  real  ocean  environments.  Modeling  provides  an 
efficient  means  by  which  to  parametrically  investigate  the  projector  (source)  radiated 
sound  pattern  under  the  influence  of  environmental  conditions  A  perfect  model  is  one  that 
perfectly  represents  reality  In  practice,  such  a  perfect  model  would  be  as  complex  as  the 
problem  it  is  attempting  to  represent  So,  there  are  some  assumptions  built  into  every 
computer  model  This  thesis  uses  a  version  of  the  PE  model  developed  at  the  University  of 
Miami  known  as  UMPE  (Smith  and  Tappert.  1994) 

In  order  to  investigate  the  effects  of  bottom  interface  roughness  and  volume 
fluctuations,  we  simulate  broadband  pulse  propagation  for  four  different  scenarios  as 
follows: 

•  flat  bottom  with  homogeneous  sediment: 

•  rough  bottom  interface  with  homogeneous  sediment; 

•  tlat  bottom  interface  with  realistic  sediment  sound  speed  fluctuations,  and 

•  rough  bottom  interface  with  realistic  sediment  sound  speed  fluctuations 


II.  BACKGROUND 

A.         THE  ACOUSTIC  MODEL 

Because  propagation  ranges  are  typically  much  larger  than  ocean  depths,  it  is 
appropriate  to  use  a  cylindrical  coordinate  system  for  long  range  propagation  Customarily 
we  define  z  as  depth,  r  as  range  outward  from  the  origin,  and  <p  as  the  azimuthal  angle 
{0—>27t)  defining  direction  along  a  radial  in  the  horizontal  plane  as  depicted  in  Figure  3  In 
general,  full  3-D  modeling  of  the  propagation  is  seldom  done  due  to  mathematical  and 
computational  complexity  It  is  often  a  good  approximation  to  the  real  ocean  acoustic 
environment  to  consider  it  azimuthally  invariant,  i  e  ,  to  assume  that  environmental 
parameters  such  as  sound-speed  profile,  water  depth,  and  bottom  composition  are 
independent  of  azimuthal  angle  to  the  source.  With  this  assumption,  we  concentrate  the 
analysis  in  depth  versus  range. 

As  mentioned  in  the  introduction,  the  PE  method  is  used  in  this  thesis  to  model 
acoustic  propagation.  The  parabolic  equation  method  was  introduced  into  underwater 
acoustics  in  the  early  1970s  by  Hardin  and  Tappert  (1973)  who  devised  an  efficient 
numerical  solution  scheme  based  on  fast  Fourier  transforms.  We  will  derive  a  PE  solution 
in  a  2-D  constant-density  fluid  environment  where  we  have  assumed  azimuthal  symmetry 
and  hence  no  dependence  on  the  ^coordinate 

We  begin  with  the  wave  equation  for  a  constant-density  medium  in  cylindrical 
coordinates  (r,z)  independent  of  <fi 

1   dlp 
c    a 

We  now  assume  we  have  a  steady  state,  time-harmonic  acoustic  field 

p(r,z,zut)  =  p(r,zy"nt  .  (2) 
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Figure  3.  Cylindrical  coordinates  in  ocean. 


Substituting  this  into  Eq.  (1)  leads  to  the  Helmholtz  equation, 

c? p      1  dp     d1 p     ,22 

^  +  ~f  +  ^T  +  kon  P  =  °  ,  (3) 

or        r  or      oz 

where    p(r,z)    is    the    acoustic    pressure,     ku  =  —     is    the    reference    wave    number. 


n{i\z)  =    ,   '   ,    is  the  index  ot  retraction,  c0is  the  reference  sound  speed  and  c(r,z)  is  the 
c\i\z) 

acoustic   sound    speed     It    is   within    c(r,z)    that    all    features   of  the   environment    are 
represented 

Tappert  (1977)  assumed  the  solution  of  Eq.  (3)  in  terms  of  an  outgoing  Hankel 
function  with  an  envelope  function  y/(r.z). 

p{r,z)  =  y/(r,z)Hi;)(k0r)  ,  (4) 

where  y/(r,z)  is  assumed  to  be  slowly  varying  in  range.  In  the  far-field    (k0r»  1),  the 
Hankel  function  can  be  approximated  by 

thus  an  alternative  relationship  between  the  acoustic  pressure,  p(i\z),   and  the  slowly 
modulating  envelope  function,  y/(i\z),  is 

p(r,z)  =  P0^y(r,z)e,k°r  .  (6) 

This  is  the  standard  definition  of  the  so-called  "PE  field  function"  y/  scaled  such  that  at 
/•  =  /<,,  |(//|  =  1  and  \p\  -  P0    Substituting  (6)  into  (3)  yields 

— -  +  2ik, ,  — -  +  — r  +  A-„  (//     -  l)y/  =  0  (7) 

a"  or       dz~ 

Now  we  introduce  the  crucial  paraxial  approximation  based  on  the  slowly  modulating 

character  of  y/, 

J-«*hf  (S) 


We  then  obtain  the  following  wave  equation, 

2tkt  — -  + — -  +  k: (tr  -  1  t//  =  0  ,  (9) 

"  cr       dz-  y  ' 

which  is  the  standard  parabolic  equation  This  can  be  rewritten  as 

V  .  (10) 


dy     ikj  ,  1    d1 


a-'  iKn  -l  +  *Z#. 

which  is  the  form  first  introduced  to  the  underwater  acoustics  community  by  Tappert 
(1974) 

Equation  (10)  is  a  one-way  wave  equation  which  can  be  solved  by  a  range- 
marching  numerical  technique.  The  UMPE  model  applied  in  this  thesis  uses  the  split-step 
Fourier  (SSF)  method  (Hardin  and  Tappert,  1973)  which  solves  the  parabolic  equation  by 
Fourier  transform  techniques  Define  the  complex  Fourier  transform  pair  as 

(//(/\r)  =  f  y/{i\k)e'k-:dk.  =  F  \v(i\k)\  ,  (II) 


M_-)  =  7~f uV{r^y'k::dz  =  l\y/{r,z)\  .  (12] 


where  k-  is  the  vertical  wave  number  We  may  rewrite  Eq.  ( 10)  as 

cty 

1) 

where 


7  =  ^o(^+^V,  d3) 


1       d1 
T    =  —7- r   ,  (14) 

op     2k;  dz2  ' 

represents  the  kinetic  energy  operator  and 

Uop=\(r>2-\)   ,  (15) 

represents  the  potential  energy  operator  This  representation  of  the  operators  as  kinetic 
and  potential  energy  operators  is  especially  insightful  when  one  wishes  to  form  the  ray 
equations  which  have  Hamiltonian  form  (Smith    1992) 

In  the  PE/SSF  method,  solutions  of  i//  are  marched  out  in  range  by 


Y(r  +  Ar)  =  e-'k^T°>e-h«ArU°>v(r)    .  (16) 

The  operator  Uop  is  simply  a  multiplication  operator  in  z-space  but  Top  contains  cross 
terms  which  couple  different  depths  However,  Top  is  a  multiplication  operator  in  vertical 
wave  number  space  It  is  therefore  desirable  to  separate  the  application  of  each  operator, 
one  in  z-space  and  one  in  /r-space  (here,  k  is  the  vertical  wave  number)  In  other  words, 
the  PE  field  function  y/  is  specified  at  some  range  r  initially  in  the  z-domain  A 
transformation  is  made  to  the  /r-domain  followed   by  a  multiplication  of  the  A--space 

operator  e  '  ''  "''  where  Top  is  the  corresponding  kinetic  energy  operator  in  the  wave 
number  domain  The  result  is  then  transformed  again  to  the  z-domain  and  is  followed  by  a 
multiplication  of  the  z-space  operator  e  '  The  final  result  is  the  field  function  at 

/•    Ai\  that  is, 

y/(r  +  A/-,z)  =  e"k^u-(r-:)  x  /r-'i*-*^'-*)  x  [f(  <//(/,  z))]}  (17) 

This  is  the  general  algorithm  of  the  PE/SSF  method. 

Previously,  we  have  assumed  the  operators  took  the  forms  defined  by  Eq.  (14)  and 
(15)  These  forms  are  commonly  referred  to  as  the  "standard  PE"  and  are  only  one  set  of  a 
number  of  various  operator  forms  In  our  simulation,  we  will  use  a  higher  order  form 
which  was  introduced  by  Thomson  and  Chapman  (1983)  referred  to  as  the  "wide-angle" 
PE  approximation  (WAPE)  with  operators  given  by 


fWAPM  =  \- 


-! 

kn  J 


(18) 


Uuir,(z)  =  -[n(z)-\]  (19) 

We  then  need  to  define  the  initial  conditions  for  the  PE  field  function,  y/(r=0,z).  It 
is  desirable  to  begin  the  calculation  by  specifying  the  source  in  the  /r-domain.  Thomson 
and  Bohun  (1988)  have  shown  that  a  wide  angle  source  can  be  defined  as 


'2jo-;(.  r*r 


^*>=^Ug  I'-bd    -(w<°-       <20) 


where  k  =k(,  corresponds  to  a  propagation  angle  of  90° 

B.  TRAVEL  TIME  ANALYSIS 

The  calculation  of  the  time  domain  arrival  structure  at  a  sjiven  range  is  an  integral 
part  of  the  analysis  in  this  thesis.  The  pressure  field  in  the  time  domain  is  represented  by 

p{r,zj)  =  F-\Po){r,z)]  =  I^M^'  .  (21) 

where  pjj'^z)  is  a  single  frequency  component.  The  arrival  time  structure  at  some  fixed 

range  r=R  can  be  computed  using  the  complex  field  p(0(R,z)  for  many  frequencies  and 

then  Fourier  transforming  to  obtain p(R, 2,/),  the  set  of  complex  pressure  values  in  the 
time  domain  To  avoid  large  transform  sizes,  the  UMPE  model  basebands  the  solutions 
and  shifts  the  center  frequency  bin  to  zero  To  reduce  the  effect  of  side-lobes,  we  apply  a 
Hanning  filter  over  the  bandwidth  before  transforming  to  the  time  domain    Additionally, 


K 

uu  — 

iknR 


we  neglect  the  phase  factor  e'  °  -  e  c"  by  simply  performing  the  Fourier  analysis  on  the 
field  function  y//0(R,z)  Therefore,  the  time  domain  is  heterodyned  around  the  value 
i,,  -R  c0  and  arrival  times  are  computed  as  values  of  "reduced  time"  or  (t-t0). 

C.         ROUGH  BOTTOM  INTERFACE 

The  ocean  bottom  roughness  is  assumed  to  be  characterized  by  a  2-D  power  law 
spectrum  at  high  wave  numbers  which  can  be  written  as 

W,\  k» !  =ak~fl  ,  (22) 

A  LcorrJ 

where  Leun  is  the  correlation  length  of  the  roughness  and  /?  is  the  spectral  exponent  (Note 
that  our  use  ot  wave  number  is  now  with  respect  to  range  scales  of  interface  roughness.) 
To  give  the  full  spectrum  a  realistically  smooth  structure,  we  assume  the  spectral  form 
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W,(k)  = —  (23) 

0  +  OV 

The  rms  roughness  represents  the  average  amplitude  of  the  undulated  bottom  interface 
variation  and  the  correlation  length  can  be  thought  of  as  a  scale  of  self-similarity  The 
normalization  factor  pi  is  defined  in  terms  of  the  rms  roughness  a,  i.e., 

27r\w2(k)kdk  =  (j2  (24) 

(i 

and  therefore 

n  ■■  2 
We  only  need  the  1-D  roughness  spectrum  along  the  slice  in  the  two  dimensional 
profiles  (depth  z  and  range  r).  Along  the  direction  of  range  propagation  in  cylindrical 
coordinates,  we  have 

W,{k)  =  l]    ,    f        W2{K)dK  ,  (26) 

k  V K~  -  k~ 


where  K  is  the  total  horizontal  wave  number  and  -k  is  the  wave  number  along  the 
propagation  direction  Substituting  Eq  (23)  and  Eq  (25)  into  Eq.  (26),  we  obtain 

W](k)  =  ya2Lcon(l  +  I:corrk2y[~;  ,  (27) 

where 


71  V  1 


\2 )    \2J  \  : 


2  f /?      ^r       dl  \2      J     \2J    \    2    J 

Hf-'JJT— ?r=~ — 7/^r  <28) 


and  /TjtJ  is  the  gamma  function 

To  avoid  the  annoying  gamma  function,  the  UMPE  model  assumes  the    1-D 
spectrum  has  the  form 

M*)=(i+02r2j  (29) 
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The  roughness  realization  is  then  simply  rescaled  by  the  desired  value  of  a  to  provide  the 
appropriate  rms  value 

D.         SEDIMENT  COMPOSITION 

Yamamoto  (1995)  summarizes  the  transformation  between  three-dimensional 
wave  number  spectra  and  the  velocitv  variabilities  which  determine  the  ocean  sediment 
properties  This  transformation  can  be  simplified  for  our  2-D  case.  The  2-D  compressional 
wave  velocity  structure  in  the  sediment  may  be  written  as 

c(x,z)  =  c0(l  +  £c(x,z))  (30) 

where  c(x,z)  is  the  velocity  at  a  point  (x,z),  c0  is  the  average  velocity  at  the  seafloor  and 
SJx.z)  is  the  relative  variability.  The  relative  variability  is  modeled  as  a  zero-mean 
quasistationary  random  function,  i.e., 

(Se(xtz))  =  0  (31) 

(Sc(x,z)Sc(x',/))  =  f\  J"  S(kxX)xcos{kx(x  -  x')  +  kz(z  -  z'))dkdk:  (32) 

The  2-D  power  spectrum  of  the  relative  velocity  variability  is  denoted  by  S(kx,k?J  By 
definition,  the  variance  [S  ])  represents  the  total  energy  of  the  power  spectrum. 

(S; )  -  f  [  s(kx ,  k.  )dkxdk.    .  (33) 

It  is  a  plausible  assumption  that  the  vertical  power  spectrum  may  be  given  by  a  power-law 
form 

St(k,)  =  Bk;fi  (34) 

where  /?  is  the  spectral  exponent  and  B  is  the  spectral  strength  constant.  The  2-D  power 

spectrum  given  by  Yamamoto  is 

dBB  i  x-(/?+|)/ 

S(k x.,k  )  =  -^—  \a-k ;  +k.:)  (35) 

4 

where  a  is  the  aspect  ratio  and  /?  2  must  be  an  integer 
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To  produce  a  random  realization  of  sediment  volume  fluctuations,  we  assume  the 
complex  wave  number  spectrum  has  a  random  amplitude  and  phase.  Ac^,  for  each  wave 
number  component  The  random  amplitude  of  each  component  of  the  power  spectrum  can 
then  be  obtained  by  A2  =-ln(r)  where  r  is  an  independent,  uniformly  distributed  random 
variable  in  the  interval  (0,1).  The  random  phase  <f>  is  another  independent,  uniformly 
distributed  random  variable  in  the  interval  (0,27t).  To  avoid  infinite  magnitudes  at  low 
wave  numbers,  the  spectrum  is  high-pass  filtered  The  filtered  spectrum  is  then  Fourier 
transformed  to  obtain  a  realization  of  the  sound  speed  variability  The  entire  field  of 
sediment  fluctuations  is  first  computed  on  a  fine  grid  and  then  interpolated  to  match  the 
computational  grid  of  the  PE/SSF  algorithm  The  entire  scheme  can  be  described  as 
follows: 


Step  1.  Construct  the  2-D  wave  number  field  [kr  (//),  k.  (/;/)) ,  according  to 

(36) 


In 
M..  = 


r 

max 


M.  = 


kr(n)  =  (n-l)Akr        n  =  l,2 N 

In 


\~max         Zbot)  (37) 

k(m)  =  (m  -  l)Mz        m  -  1,2, ...,M-  Mbot  , 
where 

'mas  =  maximum  propagation  range. 
N  -  FFT  size  in  range  direction. 
zn     =  maximum  compution  depth 
zhiil  ~  average  bottom  depth 
M  -  FFT  size  in  depth  direction, 
M .     -  index  of  average  bottom  interface 

Step  2.  Generate  two  random  numbers  series  for  the  random  phase  and 
amplitude, 

A(n,m)  =  ^-\n(anm)  ,  (38) 
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<f{n,m)  =  2nbnm  .  (39) 

where 

a,b  -  uniform  random  number  distribution  0  — >  1 . 

•  Step  3.  Compute  the  random  phase  and  amplitude  for  every  mesh  point  in  the 
2-D  realized  wave  number  field.  S'(k  {ti)k(m)) ,  i.e., 

S'(k,  (//),A-;  (//;))     s(k,(n)J:(m))A{njj,y{"",]  .  (40) 

•  Step  4.  Add  high-pass  filter  to  normalize  the  exponent  amplitude  at  low  wave 
number 

•  Step  5.  Compute  2-D  FFT  to  yield  relative  velocity  variability, 

8c{r,z)  =  2DFFT\s'{kr{n)kXm))\ . 

The  velocity  structure  of  the  sediment  realization  is  then  given  by: 

c(r,z)  =  cn+cn  xRej^/,:)).  (41) 

E.  BEAMFORMING 

In  underwater  acoustic  propagation,  much  of  the  information  about  the 
environment  and  the  relative  source/receiver  locations  is  contained  in  the  vertical  structure 
of  the  acoustic  field.  This  is  especially  true  in  shallow  water  environments  where  the 
vertical  structure  can  be  defined  in  terms  of  a  relatively  small  number  of  modes  of  the 
acoustic  waveguide  Because  these  modes  are  generally  dispersive,  the  ability  to  resolve 
the  vertical  structure  of  the  time  arrival  at  the  receiver  can  lead  to  estimations  of  the 
corresponding  source  location. 

Compared  to  the  directivity  of  single  sensors,  a  beamformer  permits  one  to  listen 
preferentially  to  wavefronts  propagating  from  one  direction  over  another.  In  this  sense,  a 
beamformer  implements  a  spatial  filter  (see,  for  example,  Defotta  et  al.,  1988,  Ziomek, 
1 995)  Consider  a  line  array  of  equally  spaced  elements  positioned  along  the  y-axis  in  an 
infinite,  homogeneous  medium  with  an  undetermined  number  of  remote  sources.  Let  the 
output  of  an  element  located  at  the  origin  of  coordinates  due  to  the  llh  source  be  Si(i). 
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Under  the  assumption  of  plane-wave  propagation  in  the  far-field,  a  source  from  direction 
0,  produces  the  following  sensor  outputs 


(42) 


where  n  is  the  array  element  index,  d  is  the  spacing,  6e  is  the  angle  of  arrival  from  the  Ith 

source,  and  c  is  the  speed  of  propagation  Beamforming  is  accomplished  by  applying 
weights  and  time  delays  to  the  individual  element  signals,  and  then  coherently  summing 
them   The  beamformer  output  when  viewed  in  the  direction    $m  is  then 

//t/(sin(#)-sin(0„,))N 


bm{i)  =  T  wnsf 


(43) 


where  the  vr„  are  weights  applied  to  each  element's  signal 

The  directional   response  characteristics  of  the  beamformer  can  be  derived  by 
considering  a  complex,  plane-wave  cw  signal  propagating  across  the  array, 

s,{t)  =  e"u<'  .  (44) 

The  source  produces  the  following  sensor  outputs  due  to  the  arrival  angle 

en(t)  =  ei<0lteink'dsm^  ,  (45) 

In      (ot 
where  ke  =  — —  -  —  is  the  wave  number  of  the  Ith  source.  The  beamformer  output  is  then 

As  C 


bnXt)  =  e'°',:Zwne'nk'^m{°f]  "M  -  W(e()e"""  , 


(46) 


where 


w(ef)  =  Yjwne 


t  mi     / 1  - .  m  I!,   i    sinfrt„p 


(47) 


represents  the  spatial  Fourier  transform  of  the  array  weights. 

A  two-dimensional  FFT  (spatial  and  temporal)  can  be  used  to  implement  the 
beamforming  operation  for  an  equally  spaced  line  array.  For  example,  consider  a  line  array 
having  N-X  equally  spaced  elements  as  depicted  in  Figure  4.  The  signals  from  each 
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element  are  sampled  simultaneously  at  rate/,  to  obtain  a  multichannel  time  series  xnm.  M- 
point  FFTs  are  taken  along  each  channel  to  obtain  the  temporal  frequency  response,  X„(j), 
where  /  is  the  frequency  bin  number,  i  e  , 

Xn(j)=TdwMx„Me'{  ,0<j<M-\  (48) 

n:    II 

The  X„(j)  are  represented  in  the  figure  as  N  horizontal  rows  of  M  FFT  bins. 

Next,  the  N  cross-channel  FFT  is  performed  with  appropriate  weighting  along  the 
array,  i.e., 

X(j,k)  = 1  »Mjy{2%)nk  ,0<k<N-\  (49) 

The  result  of  the  -V  cross-channel  FFT  series  to  spatial  frequency  (or  wave  number) 
provides  the  phase  difference  information  in  depth  for  specific  temporal  frequencies.  The 
relationships  between  the  various  indices  and  the  variables  used  in  the  beamformer  output 
are 

/,=;£/,.  (50) 

k   In 
N    d 


c   k^ 
0_.  -  sin 


(52) 


where 7  and  k  are  the  indices  in  the  frequency  and  wave  number  domains,  respectively.  In 
practice,  6m  is  determined  by  interpolating  between  values  of  ki  .  After  the  data  has  been 
interpolated  to  arrival  angle/frequency  space,  an  inverse  FFT  is  taken  along  each  angle  bin 
to  provide  arrival  angle  as  a  function  of  travel  time 
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Figure  4.  Beamforming  example. 

From  PeFotta,  1988] 
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III.  MODEL  SIMULATIONS 

As  was  described  in  the  introduction,  the  principal  characteristic  of  most  shallow- 
water  propagation  is  dominated  b\  bottom-interacting  paths  Sound  transmission  in 
shallow  water  has  an  optimum  frequency  oi  propagation  at  longer  ranges  as  shown  by 
Jensen  and  Kuperman  (1983)  Optimum  frequency  is  a  general  feature  of  ducted 
propagation  in  the  ocean.  It  occurs  as  a  result  of  competing  propagation  and  attenuation 
mechanisms  at  high  and  low  frequencies.  In  the  high-frequency  regime  there  is  increasing 
volume  and  scattering  loss  with  increasing  frequency  At  lower  frequencies  the  situation  is 
more  complicated  With  increasing  wavelength  the  efficiency  of  the  waveguide  to  confine 
sound  decreases  (the  cutoff  phenomenon)  Hence  propagation  and  attenuation 
mechanisms  outside  the  waveguide  (in  the  seabed)  become  important  In  fact,  the 
increased  penetration  of  sound  into  a  lossy  seabed  with  decreasing  frequency  causes  the 
overall  attenuation  of  waterborne  sound  to  increase  with  decreasing  frequency  Thus  we 
get  high  attenuation  at  both  high  and  low  frequencies,  while  intermediate  frequencies  have 
the  lowest  attenuation  As  Jensen's  study  shows,  the  optimum  frequency  is  in  the  range 
200-800  Hz  for  a  water  depth  of  200  m  and  we  will  use  200  Hz  as  our  center  frequency 

We  assume  that  the  ocean  bottom  sediment  is  a  fluid-like  half-space  which  means 
that  it  supports  only  one  type  of  sound  wave,  i.e.,  compressional  waves  The  information 
required  for  model  simulations  should  include  the  following  depth-dependent  material 
properties:  the  sound  speed  profile  in  the  water  column  and  compressional  wave  speed, 
attenuation  and  the  density  in  the  seabed  For  all  cases,  we  will  assume  the  water  column 
is  homogeneous  with  a  constant  sound  speed 

This  chapter  presents  the  model  simulation  results  of  four  test  cases.  Case  1 
examines  the  propagation  of  a  broadband,  impulsive  source  over  a  flat,  homogeneous 
bottom  Case  2  examines  the  propagation  in  a  similar  environment  with  a  rough  bottom 
interface  We  then  examine  which  parameters,  i.e.,  rms  roughness,  correlation  length  or 
ratio  of  them,  dominates  the  beam  resolution    Case  3  is  also  similar  to  Case  1  but  with 


realistic   bottom  volume  fluctuations  The  factors  which  most  strongly  influence  resolution 
in  Case  2  and  Case  3  will  be  combined  into  Case  4 

The  basic  range-independent  parameters  used  in  the  UMPE  model  are  listed  in 
Table  1. 


Model  Type 

Wide  Angle  PE 

Source  Function  Type 

Wide  Angle  Source 

Reference  Sound  Speed 

1500.00  m/s 

Bottom  Depth 

200.0  m 

Maximum  Computational  Depth 

400.0  m 

Mesh  Points  For  Positive  Depths 

512 

Maximum  Range 

10  0  km 

Step  In  Range 

2.0  m 

Mesh  Points  In  Range 

5000 

Source  Depth 

100.0  in 

Center  Frequency 

200.0  H/ 

Frequency  Bandwidth 

100.0  Hz 

Number  Of  Frequencies 

256 

Sound  Speed  In  Water 

1500  m/s 

Sound  Speed  In  Sediment 

1600  m/s 

Density  In  Sediment 

1500  kg/m3 

Attenuation  In  Sediment 

0.2  dB/m/kHz 

Table  1.  UMPE  Model  Parameters 

A.         CASE  1:  RANGE  INDEPENDENT 

We  begin  with  the  simple  case  consisting  of  a  constant  speed  water  column  and 
fluid-like  homogeneous  sediment  separated  by  a  flat  interface  The  corresponding 
transmission  loss  field  for  the  center  frequency  is  depicted  in  Figure  5.  We  note  the 

symmetric  constructive  interference  due  to  the  finite  number  of  trapped,    propagating 

500 


normal  modes     By  Snell's  law,  there  is  a  critical  ancle  at  arc  cos - 


600 


=20.4°  (relative  to 


horizontal)  below  which  no  energy  is  transmitted  into  the  bottom  From  the  UMPE  output 
field  with  beamforming  at  10,000  m,  we  obtain  the  angle-time  plot  shown  in  Figure  6  As 
expected,  this  energy  arrives  primarily  near  horizontal  with  an  aperture  of  approximately 
±20°  The  steeper-angle  energy  >20°  is  rapidly  lost  by  transmission  into  the  seabed. 
Furthermore,  the  lowest  modes  arrive  first,  before  the  higher  modes,  as  would  be  expected 
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Figure  5.  Range-depth  transmission  loss  field  for  flat  bottom  interface. 
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Figure  6.  Arrival  angle-time  plot  for  flat  bottom  interface. 


from  an  analysis  of  the  modal  dispersion  in  this  environment  A  standard  calculation  of  the 
number  of  trapped  modes  was  found  to  be  consistent  with  the  number  of  arriving  modal 
peaks  when  the  range  of  the  transmission  loss  palette  was  extended  to  expose  the  weakest 
modal  arrivals. 

B.         CASE  2:  ROUGH  BOTTOM  INTERFACE 

A  rough  seafloor  will  cause  attenuation  of  the  mean  acoustic  field  propagating  in 
the  ocean  waveguide  due  to  incoherent  scattering  thereby  decreasing  the  mean  signal  to 
noise  ratio  For  individual  pulse  transmissions,  the  instantaneous  signal  becomes  smeared 
in  time  and  space  due  to  mode  coupling  A  rough  bottom  interface  is  added  to  the 
environment  of  Case  1  in  order  to  investigate  the  effect  of  scattering  on  the  instantaneous 
resolution  of  the  pulse  in  the  water  column 

As  described  in  Chapter  11,  there  are  two  parameters  that  control  the  roughness, 
i  e  ,  the  nils  roughness  and  the  correlation  length  In  order  to  determine  the  dominant 
factor,  we  simulated  two  sets  of  comparisons  In  one,  the  correlation  length  is  fixed  at 
1000  m  and  the  rms  roughness  has  values  of  2,  5,  and  10  m  In  the  other,  the  rms 
roughness  is  fixed  at  2  m  with  correlation  lengths  of  40,  200,  and  400  m  Figure  7-12 
show  the  beamformed  angle-time  plots  at  10,000  m.  We  can  see  that  increasing  the  rms 
roughness  and  decreasing  the  correlation  length  will  reduce  the  resolution  but  it  is 
difficult  to  see  which  factor  dominates  However,  the  distortion  changes  look  similar  and 
suggest  that  the  ratio  of  rms  roughness  and  correlation  length  may  be  the  important 
parameter  We  simulated  several  realizations  with  fixed  ratio  of  the  rms  roughness  and  the 
correlation  length  to  confirm  our  hypothesis,  specifically,  (rms  roughness):  (correlation)  = 
1:500,  1  100  and  1:20  The  beamformed  results  are  displayed  in  Figure  13-15.  Combined 
with  the  previous  simulations,  we  find  that  higher  ratios  retain  higher  resolution  and  lower 
ratios  degrade  the  resolution  more  severely  Figure  16-18  show  single  frequency  (200  Hz) 
transmission  loss  fields  in  depth-range  slices  for  various  rough  interface  realizations. 
(  ompared  to  the  flat  bottom  case,  we  see  that  there  is  more  energy  penetration  into  the 
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Figure  7.  Arrival  angle-time  plot  for  rough  bottom  interface  with  rms  roughness  = 
2m  and  correlation  length  =  1,000m. 
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Figure  8.  Arrival  angle-time  plot  for  rough  bottom  interface  with  rms  roughness  = 
5m  and  correlation  length  =  1,000m. 
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Figure  9.  Arrival  angle-time  plot  for  rough  bottom  interface  with  rms  roughness 
10m  and  correlation  length  =  1,000m. 
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Figure  10.  Arrival  angle-time  plot  for  rough  bottom  interface  with  rms  roughness  = 
2m  and  correlation  length  =  40m. 
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Figure  11.  Arrival  angle-time  plot  for  rough  bottom  interface  with  rms  roughness 
2m  and  correlation  length  =  200m. 


27 


s 

o 

o 

<* 

II 

f 

a 

i—i 

c 

O 

a 

1 

s 

<M 
II 

to 


CM 

O 


u 

to 
a) 

s 


13 

a> 
o 


o 

GO 


CD 


Tt< 


<N 


pq 


(Sap)  9{§UB  {BALLiy 


Figure  12.  Arrival  angle-time  plot  for  rough  bottom  interface  with  rms  roughness 
2m  and  correlation  length  =  400m. 
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Figure  13.  Arrival  angle-time  plot  for  rough  bottom  interface  with  rms  roughness  = 
5m  and  correlation  length  =  2,500m. 
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Figure  14.  Arrival  angle-time  plot  for  rough  bottom  interface  with  rms  roughness 
2m  and  correlation  length  =  200m. 
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Figure  15.  Arrival  angle-time  plot  for  rough  bottom  interface  with  rms  roughness 
5m  and  correlation  length  =  100m. 


31 


L  O 


E 
o 
o 

LO 

c\f 

ii 

■ 

o 
o 

E 

to 

II 

C/) 

E 

CM 

0 
CO 
03 

o 


E 
o 

C 
CO 


o 


(ai)  md9Q 


Figure  16.  Range-depth  transmission  loss  field  for  rough  bottom  interface  with  rms 
roughness  =  5m  and  correlation  length  =  2,500m. 
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Figure  17.  Range-depth  transmission  loss  field  for  rough  bottom  interface  with  rms 
roughness  =  2m  and  correlation  length  =  200m. 
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Figure  18.  Range-depth  transmission  loss  field  for  rough  bottom  interface  with  rms 
roughness  =  5m  and  correlation  length  =  100m. 


bottom  at  longer  ranges,  and  the  stable  interference  structure  has  broken  down  into  a 
complex  pattern. 

C.         CASE  3:  SEDIMENT  VOLUME  FLUCTUATIONS 

The  scattering  due  to  the  velocity  variability  within  the  sediment  is  considered  to 
be  the  major  scattering  mechanism  in  many  geological  fields  in  shallow  water  based  on 
extensive  measurements  of  acoustic  backscattering  and  sediment  cores  (Jackson  and 
Briggs  1992,  Yamamoto  1995).  The  geology  of  shallow-water  seabeds  is  complicated 
because  of  the  physical,  chemical,  and  biological  activities  which  took  place  though  the 
long  geological  time  spans.  The  measurements  of  the  velocity  and  density  variabilities  in 
the  ocean  are  beyond  this  thesis  subject.  For  our  simulations,  we  shall  use  the  sediment 
model  described  in  Chapter  II  with  some  typical  values  extracted  from  Yamamoto's 
work. 

In  this  case,  we  use  a  flat  bottom  interface  and  add  the  sediment  volume 
perturbation.  According  to  Yamamoto's  model,  the  bottom  sediment  perturbation  is 
controlled  by  three  parameters:  spectral  exponent  P,  aspect  ratio  a,  and  the  spectral 
strength  constant  B.  In  his  conclusions,  the  spectral  exponent  P  and  aspect  ratio  a  are 
fairly  constant  regardless  of  the  sediment  type  but  the  spectral  strength  constant  B  varies 
depending  on  the  sediment  type  and  location.  We  will  fix  a  and  P  and  vary  the  values  of 
B  by  choosing  three  typical  spectral  strength  constants,  e.g.,  2x10'  ,  2x10"  and  2x10" 
(units  of  B  are  m  "  ).  The  sediment  field  for  strength  constant  B  =  2x10"  m  "  is  depicted 
in  Figure  19.  The  rms  value  of  the  sound  speed  fluctuations  in  the  sediment  for  this  value 
of  B  is  22.4  m/s.  The  beamformed  results  for  the  three  values  of  B  are  shown  in  Figure 
20-22.  The  pulse  shape  does  appear  to  degrade  slightly  as  the  spectral  strength  constant  is 
increased.  Overall,  however,  the  sediment  volume  perturbation  doesn't  seem  to  affect  the 
pulse  resolution  significantly. 
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Figure  19.  Sediment  field  for  strength  constant  B=2xl0    in 
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Figure  20.  Arrival  angle-time  plot  for  sediment  volume  fluctuation  with  strength 


constant  B=2xl0'3m1"n. 
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Figure  21.  Arrival  angle-time  plot  for  sediment  volume  fluctuation  with  strength 
constant  B=2xl0"4 mlp. 
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Figure  22.  Arrival  angle-time  plot  for  sediment  volume  fluctuation  with  strength 
constant  B=2xl05mlp. 
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D.  CASE  4:  ROUGH  BOTTOM  INTERFACE  WITH  SEDIMENT  VOLUME 
FLUCTUATIONS 

We  have  seen  the  effect  of  rough  bottom  interfaces  and  sediment  volume 
fluctuations  separately  We  now  need  to  combine  these  two  factors  From  Case  2  and 
Case  3  results,  it  appears  that  the  effect  of  a  rough  bottom  interface  dominates  the 
pulseresolution  The  strongest  sediment  volume  perturbation  strength  constant,  B  =  2x  10° 
/w  ,  was  chosen  to  produce  the  largest  volume  effect.  The  interface  correlation  length 
was  fixed  at  1000  m  and  the  rms  roughness  was  varied  from  2  m  to  10  m.  Comparing 
these  results  to  those  cases  with  no  volume  fluctuations  showed  little  difference  indicating 
that  these  levels  of  interface  roughness  dominate  the  effects  on  pulse  resolution. 

E.  PULSE  RESOLUTION  ANALYSIS 

We  have  shown  that  the  effect  of  the  bottom  roughness  and  sediment  volume 
fluctuation  degrades  low  frequency  pulse  resolution  in  shallow  water  propagation.  In  this 
section,  we  will  quantify  the  effect  of  rough  bottom  and  sediment  volume  fluctuation  on 
pulse  resolution  in  order  to  obtain  a  clearer  relationship.  The  conventional  angle  resolution 
is  defined  as  an  angle  width  of  Me  peak  amplitude  In  order  to  simplify  the  Case  2  analysis, 
we  fixed  the  correlation  length  at  a  value  of  1000  m  and  varied  the  rms  roughness  Figure 
23  shows  several  slices  in  reduced  times  of  I,  15,  2,  2.5  and  3  seconds  from  an  arrival 
angle-time  plot  for  rms  roughness  of  0  (the  Case  1  results),  2,  5,  10  and  20  m.  Each 
reduced  time  data  point  is  actually  an  average  of  five  points  spread  over  a  time  bin  width 
of  0.2  sec,  i.e.,  the  datum  for  the  reduced  time  of  1  sec  is  an  average  of  data  at  0.90,  0.95, 
1.00,  1  05,  and  1.10  sec.  Polynomial  fits  to  the  data  produce  the  smooth  curves  shown  in 
the  figures  We  can  see  a  good  angle  resolution  at  lower  rms  roughness  values  As  the  rms 
roughness  is  increased,  the  resolution  is  found  to  degrade,  particularly  for  the  later  arrival 
times  Fins  is  expected  since  later  arrivals  correspond  to  higher  modes  which  interact  more 
strongly  with  the  bottom  As  an  estimation  of  this  effect,  we  note  that  the  angular  pulse 
width  appears  to  double  for  an  rms  value  of  ~I0  in  corresponding  to  an  rms:  correlation 
ratio  of  I   100 
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In  Figure  24.  slices  are  taken  at  arrival  angles  0,  5,  10,  and  20  deg  In  this  case, 
three  points  over  a  width  of  1  deg  are  averaged  to  produce  each  data  point  shown.  There 
is  an  anomalous  point  at  an  rms  of  2  m  due  to  unresolved  peaks  of  two  modes.  This  point 
was  not  used  in  the  curve  fit.  The  time  resolution  degrades  when  the  rms  roughness 
increases  for  all  arrival  angles  except  0  deg.  This  is  due  to  less  bottom  interactions  of  the 
horizontal  propagation  Furthermore,  we  notice  that  the  time  width  does  not  degrade  as 
rapidly  as  the  angular  width  None  of  the  rms  roughness  values  appear  to  double  the  pulse 
width  in  time  1  he  peak  minimum  transmission  loss  also  increases  with  rms  roughness  as 
shown  in  Figure  25  by  roughly  5-10  dB  at  an  rms  value  of  10  m 

To  perform  a  similar  analysis  of  the  Case  3  environment,  the  spectral  exponent  and 
aspect  ratio  were  fixed.  Obviously,  the  sediment  volume  fluctuations  show  more  stable 
results  with  varying  strength  constant  as  depicted  in  Figure  26-28.  The  angle  and  time 
resolution  change  little  with  increasing  strength  constant.  The  peak  (minimum) 
transmission  loss  also  remains  stable. 

In  Figure  29,  the  peak  minimum  transmission  loss  computed  with  the  combined 
effects  of  a  rough  bottom  interface  and  sediment  volume  fluctuations  shows  very  similar 
results  as  in  Figure  25  It  is  reasonable  to  conclude  that  these  levels  of  bottom  roughness 
dominate  pulse  resolution  effects  in  simple  shallow  water  propagation 
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Figure  23.  Angle  width  versus  rms  roughness  for  reduced  time  =  1,  1.5,  2.  2.5  and  3 
second. 
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Figure  24.  Time  width  versus  rms  roughness  for  arrival  angle  =  0,  5,  10,  15  and  20 
degree. 
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Figure  25.  Peak  minimum  transmission  loss  versus  rms  roughness  for  reduced  time 
=  I,  1.5,  2.  2.5  and  3  second. 
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Figure  26.  Angle  width  versus  strength  constant  for  reduced  time  =  1,  1.5,  2,  2.5  and 
3  second. 
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Figure  27.  Time  width  versus  strength  constant  for  arrival  angle  =  0,  5,  10,  15  and 
21)  degree. 
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Figure  28.  Peak  minimum  transmission  loss  versus  strength  constant  for  reduced 
time  =  1,  1.5,  2,  2.5  and  3  second. 
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I  igure  29.  Peak  minimum  transmission  loss  versus  rms  roughness  for  reduced  time 
'-  I.  1.5,  2,  2.5  and  3  second  combine  with  strength  constant  =  2xl0"3mlp  . 
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IV.   CONCLUSION 

In  this  thesis,  we  have  presented  a  numerical  investigation  of  the  effects  of  bottom 
interface  roughness  and  sediment  volume  fluctuations  on  the  resolution  of  broadband 
pulse  propagation  in  a  simple,  shallow  water  environment.  Predictions  of  the  acoustic 
propagation  were  accomplished  with  a  broadband,  wide-angle  parabolic  equation  model. 
Fluctuations  of  the  bottom  interface  and  volume  were  computed  from  realistic  spectra 
characterizing  each  structure.  The  strength,  or  level  of  variability,  of  each  was  varied 
independently  in  order  to  determine  the  relative  influence  on  pulse  resolution 

To  be  consistent  with  measurement  techniques  employed  in  underwater  acoustics 
experiments,  a  beamforming  algorithm  was  written  to  discriminate  vertical  arrival  angles 
as  a  function  of  arrival  time  When  no  fluctuations  were  introduced,  the  solution  to  the 
range-independent  problem  was  found  to  agree  with  a  standard  normal  mode  analysis 
When  fluctuations  were  introduced,  the  symmetrical  arrival  time  pattern  became  distorted, 
indicating  a  degradation  in  the  resolution  of  the  pulse  Such  degradations  are  expected  to 
limit  the  ability  of  some  operational  naval  systems 

When  interface  roughness  was  introduced,  the  results  suggested  that  the  main 
parameter  influencing  resolution  was  the  ratio  of  rms  roughness  to  correlation  length  In 
other  words,  a  high  value  of  rms  roughness  does  not  independently  imply  a  degradation  in 
the  pulse  shape  if  the  correlation  length  is  quite  large.  Similarly,  a  relatively  small  rms 
roughness  may  still  significantly  affect  pulse  resolution  if  accompanied  by  a  very  small 
correlation  length  In  contrast,  the  largest  realistic  scales  of  volume  fluctuations  had  little 
impact  on  the  pulse  shape  A  combination  of  large  volume  fluctuations  and  various  scales 
of  rms  interface  roughness  showed  that  degradations  of  pulse  resolution  are  dominated  by 
the  interface  structure  While  it  may  be  possible  that  a  different  environment  (eg  ,  strongly 
downward  refracting  water  column  and  upward  refracting  sediment)  would  increase  the 
influence  of  sediment  heterogeneity,  such  cases  would  also  tend  to  increase  interactions 
with  the  bottom  interface  resulting  in  similar  conclusions. 
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An  attempt  was  made  to  quantity  these  results  by  measuring  the  Me  -  width  of  the 
pulse  at  various  arrival  angles  or  arrival  times  as  a  function  of  the  fluctuation  strength.  It 
was  observed  that  an  rms  interface  roughness  of  1 0  m  with  a  correlation  length  of  1 000  m 
produced  a  doubling  of  the  angular  width  of  the  arrivals  and  increased  the  transmission 
loss  by  as  much  as  10  dB  The  time  width  was  more  stable,  however  In  none  of  the  cases 
with  volume  fluctuations  was  the  signal  found  to  degrade  appreciably  When  the  two 
influences  were  combined,  the  results  were  nearly  identical  to  the  case  with  no  volume 
fluctuations 

Overall,  these  results  suggest  that  bottom  interactions  in  shallow  water 
propagation  are  dominated  by  interface  scatter  The  degradation  of  pulse  resolution 
depends  on  a  combination  of  the  rms  interface  roughness  and  an  associated  correlation 
length  scale  In  order  to  know  a  priori  the  limitations  of  pulse  resolution,  it  is  necessary  to 
obtain  an  accurate  statistical  representation  of  the  sea  floor  topography.  In  contrast,  the 
sub-bottom  volume  fluctuations  do  not  appear  to  significantly  influence  the  forward 
propagation 
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